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ABSTRACT 



The Saturnian coorbital satellites Janus and Epimetheus present a unique dy- 
namical configuration in the Solar System, because of high-amplitude horseshoe 
orbits, due to a mass ratio of order unity. As a consequence, they swap their 
orbits every 4 years, w hile their orbital periods is about 0.695 days. Recently, 



Tiscareno et al.l (120091 ) got observational informations on the shapes and the ro- 



tational states of these satellites. In particular, they detected an offset in the 
expected equilibrium position of Janus, and a large libration of Epimetheus. 

We here propose to give a 3-dimensional theory of the rotation of these satel- 
lites in using these observed data, and to compare it to the observed rotations. 
We consider the two satellites as triaxial rigid bodies, and we perform numerical 
integrations of the system in assuming the free librations as damped. 

The periods of the three free librations we get, associated with the 3 dimen- 
sions, are respectively 1.267, 2.179 and 2.098 days for Janus, and 0.747, 1.804 
and 5.542 days for Epimetheus. The proximity of 0.747 days to the orbital period 
causes a high sensitivity of the librations of Epimetheus to the moments of iner- 
tia. Our theory explains the amplitude of the librations of Janus and the error 
bars of the librations of Epimetheus, but not an observed offset in the orientation 
of Janus. 



Subject headings: Saturn, satellites — Resonances, spin-orbit — Rotational dy- 
namics 
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Introduction 



The Saturnian system has many coorbital satelhtes, hke Tethys-Telesto-Calypso and 
Dione-Helen-Polydeuces, but the pair Janus-E pimetheus presents a unique configuration be- 
cause their not-so-big ma ss ratio, i.e. ~ 3.604 (|Jacobson et al.1120081) . induces large hbrations 



named horseshoe orbits (iDermott fc Murravl Il98lh . These inner satelhtes have originally 
been observed by Dollfus in 1966 JPoUfusl l Il967l ). who thought to have discovered only one 
satellite, named Janus, but found it di fficult, if not imposs i ble, t o fit a reliable orbit of this 
body to observations. Ten years later, [Fountain fc LarsonI (119771 ) showed that this discrep- 
ancy could be explained by the presence of at least one another satellite, now known as 
Epimetheus, whose exist ence was confi r med i n 1980 thanks to Earth-based observations and 
Voyager images (see e.g. iLarson et al.l (Il98ll )). 



The horseshoe orbits are the consequence of a 1 : 1 orbital resonance (i.e., Janus and 
Epimetheus have the same mean orbital period), that provokes an orbital swap every 4 
years, the two satellite recovering their original semimajor axis after a second swap. From 
a dynamical point a view, we can say that the period of libration associated with this 1 : 1 
resonance is 8 years. 

The Cassini spacecraft, currently orbiting in the Saturnian system, presents a unique 
opportunity to observe thes e orbital swaps and thei r consequences, for instance on the density 
waves in the rings (see e.g. iTiscareno et al.l (120061 )). Recently, Cassini images of Janus and 
Epimetheus allowed to derive the shapes of these bodies, an estimation of their moments of 
inertia, and also some measures of their rotation. In particular, it has been established that 
there is a permanent offset of several degrees between Janus minimum moment of inertia and 
the equilibrium sub-S aturn point, and that Epimetheus has an oscillation about synchronous 
rotation of 5.9°±1.2° (ITiscareno et al.ll2009l ). These last two results are expected to give clues 
on the distribution of masses inside of these two bodies, as did for instance the detection of 



an oscillation of Mercury around the 3 : 2 orbital resonance (IMargot et al.ll2007l ). considered 
as the evidence of a molten core. 

The goal of this paper is to compute a theory of the rotation of Janus and Epimetheus, 
and to compare it with the observational data. We first review some aspects of the orbital 
dynamics of Janus and Epimetheus, because they have a direct influence on the rotation 
through the Saturnian perturbation. In particular, we express the fundamental frequencies 
of the orbital perturbations, and interpolate JPL/HORIZONS ephemerides before including 
them into a 3-degrees of freedom numerical model of the rotation of Janus and Epimetheus, 
seen as rigid bodies. Finally, we propose some interpretations of the observed data, in 
comparison with our theory. 
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2. An analysis of the orbital motion 



It is necessary to know as accurately as possible the orbital motion of Janus and 
Epimetheus about Saturn, to estimate the influence of the Saturnian torque on their rotation. 
The most reliabl e ephemerides of Janus and Epimetheus available are JPL HORIZONS' ones 
0, elaborated by Jacobson et al. ( 2008 ). This server provides data tables of the positions, 
velocities, orbital elements,. . . of the most important Solar System objects (including Janus 
and Epimetheus) at given dates. We explain in this part how we extract from JPL data the 
useful informations that will help us to understand the orbital dynamics of these bodies. 



2.1. The proper modes 

The flrst information we are interested in is the proper modes of the orbital motion, 
i.e. the fundamental frequencies of the system. If we assume that the orbital motion of 
a given body is not chaotic, its variables can be expressed under a synthetic form, i.e. an 
inflnite sum of sinusoidal terms, that depend on a limited number of fundamental frequencies, 
this number depending itself on the number of involved bodies. If we neglect the periodic 
influence of the other Saturnian satellites and consider that our system is only composed 
of an oblate Saturn, Janus and Epimetheus, we can restrict our fundamental frequencies to 
6 elements, i.e. 3 for each satellites, these frequencies being the orbital one, and the slow 
frequencies of precessions of nodes and pericentres. The 1:1 orbital resonance forces the 
equality of the two mean motions, but also results in the apparition of a proper frequency 
associated with this resonance, of period close to 8 years, i.e. two orbital swaps. So we have 
6 proper modes. 

From the JPL cartesian coordinates (NAIF kernels SAT299 for Janus and Epimetheus, 
and SAT317 for Saturn) and masses of Janus and Epimetheus, we derived the following 
keplerian elements, in the reference frame centered on the center of mass of Saturn, and 
referring to the equatorial plane of Saturn and the node of this plane with the ecliptic at 
J2000.0: 

• a: semi-major axis 

• A: mean longitude 



^http:/ /ssd.jpl.nasa.gov/?horizons 
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z = k + \/^^h = eexp {\^—ltu), e and w being respectively the eccentricity and the 
longitude of the pericentre 



C = Q + y/~^P = sin exp (a/— 1^), / and ^ being respectively the inclination and 
the longitude of the ascending node. 



These elements are the classical elliptic ones. Another possibility could be to use the 
epicyclic elements. The basic idea is to consider the oblateness of Saturn in the computation 
of the osculating elements. The third Kepler law is being modified this way: 



3 



R 



0(J4 



(1) 



n being the instantaneous mean motion of the considered body (Jan u s or E pimetheus), 
M-^ Saturn's mass, and R-^ its equatorial radius (see e.g. iGreenbergI (jl98ll )). This for- 
mulation does not change the fundamental frequencies, because they physically correspond 
to revolutions of the body about its parent planet, but change the semimajor axes, ec- 
centricities and inclinations. In particular, it can be shown that it drastically reduces 
the amplitudes of th e short period oscillations in the case of a small body orbiting close 
to an oblate planet f Brouwerl 19591). Such a formulation is very convenient for describing 



the dynamics of planetarv rings (IBorderies fc Longarettilll987l:lLongaretti fc Borderieslll991 



Borderies-Rappaport &: Longaretti 1994 ; Renner Sz Sicardy 20061 ) because it gives quite con- 



stant elements and a smaller eccentricity than the classical elliptic elements, but in our case 
we want to identify the mean orbital period. Moreover, no influence is expected on the 
representation of the orbital swap. This is the reason why we did not use them. 

Once these elements have been obtained as data tables, we performed frequency analyses 
to get periodic time series. The frequency analysis algorithm we used is based on Laskar's 
original idea, na r ned N AFF as Numerical A nalysis o f the F undamental Frequencies (see for 
instance iLaskarl (119931 ) for the method, and iLaskarl (120031 ) for the convergence proofs). It 



aims at identifying the coefficients and of a complex signal f{t) obtained numerically 
over a finite time span [— T; T] and verifying 



/W ~ J^afcexp(v^t^fct), 



(2) 



k=l 



where LOk are real frequencies and complex coefficients. If the signal f{t) is real, its 
frequency spectrum is symmetric and the complex amplitudes associated with the frequencies 
LJk and —Uk are complex conjugates. The frequencies and amplitudes associated are found 
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with an iterative scheme. To determine the first frequency cji, one searches for the maximum 
of the amphtude of 

<f){uj)=<f{t),exp{V^ut)>, (3) 
where the scalar product < f{t),g(t) > is defined by 

1 



2T 



<f(t),git)>=— / fit)git)xit)dt, (4) 



T 



and where x{t) is a weight function, i.e. a positive function with 



±, I ^Xit)dt = 1. (5) 



Once the first periodic term exp(v^co'it) is found, its complex amphtude ai is obtained 
by orthogonal projection, and the process is started again on the remainder fi(t) = fit) — 
aiexp{^/^ujit). The algorithm stops when two detected frequencies are too close to each 
other, what alters their determinations, or when the number of detected terms reaches a 
maximum set by the user. This algorithm is very efficient, except when two frequencies are 
too close to each other. In that case, the algorithm is not confident in its accuracy and stops. 
When the difference between two frequencies is larger than twice the frequency associated 
with the length of the total time interval, the determination of each fundamen tal frequency 



i s not perturbed by the other ones. Although the iterative method suggested by lChampenois 



( 1l998l ) allows to reduce this distance, some troubles still remain when the frequencies are 



too close to each other. 

We used ephemerides over the time interval 1950-2050, i.e. the widest interval over 
which the JPL ephemerides are available. This time span is quite long in comparison with 
the orbital period of these two satellites (~ 17 hours), and long enough for estimating the 
period of the resonant argument 8 years, i.e. two orbital swaps). The shortest unaliased 
period that can be detected is twice the time step. Here our time step was 3 hours, what is 
short enough to detect 8h-contributions, i.e. half the orbital period expected. 

The frequency analyses of the keplerian elements of Janus and Epimetheus allowed us 
to identify 7 proper modes (see TabU]), written as A, 0, wj, we, ^j, aiid u. The first six 
proper modes are the ones expected, while the last one, u, has been actually detected in the 
frequency analyses. This seventh proper mode is difficult to identify, particularly because its 
period is larger than our time interval, so we lack of accuracy on its period. Moreover, the 
amplitude associated is small. A is close to the mean longitude of Janus, is the resonant 
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argument, wj and we are close to the longitudes of the pericentres of Janus and Epimetheus, 
while and are close to the longitudes of their ascending nodes. By "close to", we 
mean that the proper mode is the main component of the keplerian element associated with. 



Table 1: Proper modes of the orbital motions of Janus and Epimetheus. The column " Origin" 
gives the variable from which the given numerical values of the proper mode hav e been 
extracted, while the last column gives the frequencies given by lJacobson et al.l (120081 ). fitted 
on [2003; 2005]. 





Frequency (rad/y) Phase at J2000 


Period 


Origin 


JSP2008 


A 


3304.0143278 


-114.564° 


0.69459 d 


Aj 


3303.6716 










Xe 


3305.2554 





0.7847244 


177.674° 


8.00687 y 


aj 




Wj 


13.0908741 


129.064° 


175.30788 d 


Zj 


13.0869 


We 


13.0928523 


-121.751° 


175.28140 d 


ze 


13.1022 




-13.0386776 


114.152° 


-176.00968 d 





-13.0359 




-13.0400438 


152.811° 


-175.99124 d 


Ce 


-13.0512 


UJ 


0.0461439 


-120.692° 


136.16498 y 








The numerical values we give in this table come from outputs of the frequency analyses. 
We can see that some frequencies are very close to each other, for instance for the precessions 
of the pericentres. They could be distinguished from each others in the same signal if the 
interval of study were longer than twice the period associated with the difference of their 
frequencies, i.e. ~ 200,000 years. Nevertheless, we can also use the phases to identify 
the contribution of each proper mode, that is the reason why we are confident in their 
identification. Our n umerical frequencies are in good agreement with the ones given by 



Jacobson et al.l (120081 ) . that have been obtained from a fit over two years. 



We here give the frequency analyses of the two semi-major axes (TabE] for Janus and 
Tab J3] for Epimetheus). The decomposition of the other orbital variables can be found in 
the supplemental material. These tables should be used this way: the TabjS] means that the 
semimajor axis of Janus could be approximated by 



a{t) ^ 152043.049 

+ 13.182 cos(0(t)) 

+ 7.397 cos(A(t) - roj(t)) 

- 4.538 cos(30(t)) 

+ ... (6) 
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with</)(t) = 0.7847244t+177.674°, A(t) = 3304.0143278t-114.564° and wj(t) = 13.0908741t+ 
129.064° (Tab{T]), the amplitudes being in km, the frequencies in rad/y, and the time in years, 
the origin being J2000. The cosines should be replaced by sines or by complex exponentials 
when it is stated in the caption. The cosines and sines are suitable for real variables, and 
the opportunity to use sine or cosine is decided in reading the phases. 

Table 2: Semi-major axis of Janus. The series are in cosine. 



N 


A 





Wj 


Amplitude (km) 


Period 


1 


_ 




_ 


152043.049 


oo 


2 


_ 


1 


_ 


13.182 


8.00687 y 


3 


1 


_ 


-1 


7.397 


0.69735 d 


4 


_ 


3 


_ 


-4.538 


2.66892 y 


5 


_ 


5 


_ 


2.403 


1.60135 y 


6 


1 


-1 


-1 


2.098 


0.69752 d 


7 


1 


1 


-1 


2.097 


0.69719 d 


8 




7 




-1.547 


1.14382 y 


9 




9 




1.087 


324.94105 d 


10 




11 




-0.802 


265.86068 d 


11 




13 




0.611 


224.95882 d 


12 




15 




-0.476 


194.96413 d 


13 




17 




0.377 


172.02726 d 


14 


1 


2 


-1 


0.365 


0.69702 d 


15 


1 


-2 


-1 


0.359 


0.69769 d 


16 




19 




-0.302 


153.91905 d 


17 




21 




0.245 


139.26011 d 


18 


1 


-3 


-1 


0.203 


0.69785 d 


19 




23 




-0.201 


127.15028 d 


20 


1 


3 


-1 


0.200 


0.69685 d 


21 




25 




0.165 


116.97818 d 



A striking point in these decompositions are the odd harmonics of 0. Their presence 
can be explained in considering that the semimajor axes are close to square wave, and that 
the Fourier series of a square wave is composed only of odd harmonics. We can also notice 
the quite slow decrease of the amplitude, showing that the quasi-periodic decomposition of 
such a signal converges slowly. Here, all the terms detected by the frequency analysis are 
given, there is no cut-off based on the amplitude. The estimation of the accuracy of this 
representation will be given later on the positions of the bodies (Figl2]). 
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Table 3: Semi-major axis of Epimetheus. The series are in cosine. 



N 


A 




ZD 77 


Amplitude (km) 


Period 


1 








152043 602 


oo 


2 




1 




-47.500 


8 00675 V 


3 




3 




16.353 


2 66892 V 


4 




5 




—8.660 


1 60135 V 


5 


1 


-1 


-1 


-6.251 


69752 d 


6 


1 


1 


-1 


-6.234 


69719 d 


7 




7 




5.577 


1 14382 V 


8 


1 


2 


-1 


-4.731 


69702 d 


9 


1 


-2 


-1 


-4.726 


69735 d 


10 




9 




-3.918 


324.94138 d 


11 


1 




-1 


-3.114 


0.69735 d 


12 




11 




2.889 


265.86079 d 


13 




13 




-2.201 


224.95893 d 


14 




15 




1.715 


194.96444 d 


15 




17 




-1.359 


172.02730 d 


15 


1 


3 


-1 


-1.319 


0.69685 d 


16 


1 


-3 


-1 


-1.293 


0.69785 d 


17 




19 




1.090 


153.91905 d 



In addition to these frequency analyses, we also find drifts in the orbital elements of 
Janus and Epimetheus, indicating the presence of long-term effects, that a representation 
over one century cannot render. We made a linear least-squares fit of the eccentricities e and 
inclinations I of Janus and Epimetheus and we get, for Janus, e = (3.93785 x 10~^ ± 1.708 x 
10"^)t + 7.30539 X 10-3±4.93 x 10-^ and / = ((-9.31642 x 10-^±5.263 x 10-6)t + 9.86587± 
1.519 X 10"^) arcmin, and for Epimetheus e = (-1.03189 x 10"^ ± 1.74 x 10^^)t + 1.01621 x 
10-2 ± 5.023 X 10-6 and / = ((1.49189 x lO^^ ± 5.529 x m-^)t + 21.1834 ± 1.596 x 10-^) 
arcmin, the time unit being the year and its origin J2000.0. 



2.2. Focus on the orbital swap 



The orbital swap occurring every 4 years is the key point of the dynamics of Janus and 
Epimetheus. We here p ropose to recall its main aspects, that have been ex t ensively studied 



i n pre vious works, e.g. lYoder et al.l (Il983l ): lYoder. Synnott fc Said (119891 ): [Nicholson et al. 
Jl992h . 
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1995 2000 

Date fv) 




(b) 




(c) (d) 

Fig. 1. — The orbital swap, plotted with JPL HORIZONS data. The panels (a) and (b) 
show the variations of the semi-major axis of Epimetheus, while the panel (c) shows the one 
of Janus. The panel (d) represents the variations of the mean longitude of Epimetheus after 
removal of a slope of 3304.01449932 rad/y. 



The Figure [T] shows the variations of the semi-major axes and the mean longitudes of 
Epimetheus because of the orbital swap. These figures have been obtained with the keplerian 
elements derived from the JPL HORIZONS cartesian coordinates. The most striking is the 
exact synchronization between the orbital swaps of the two bodies (panels (a) and (c)), 
illustrating the energy exchanges due to the 1 : 1 orbital resonance. We can also notice the 
amplitudes of the swaps, respectively ~ 13 km for Janus and ~ 47 km for Epimetheus. The 
ratio between these two amphtud es is the mass ratio o f the two bodies (i.e. ~ 3.6) , as already 
shown by previous authors (e.g. iDermott fc Murrayl (Il98ll ). lYoder et al.l (Il983l )). We also 
notice the thickness of the plateau, due to short period oscillations whose amplitude partly 
depends on the orbital eccentricity of the considered body. As already said, this plateau 
would have been thinner with epicyclic elements. The consequences of the resonance on the 
mean longitudes are 8-y periodic sawtooth waves. 
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Finally, the panel (b) is a zoom on the transition, similar to the Figure 2b of lNicholson et al. 



(119921 ) ■ It illustrates the motion of Epimetheus during the swap. We can graphically evaluate 
the duration of the swap at about 6 months. 



2.3. Using the ephemerides in the numerical computations 

The quasi-periodic decomposition of the orbital ephemerides is a good way to get a 
description of the ephemerides that is physically understandable, but it can be reliably 
used in numerical computation only if its residuals with the real ephemerides are small 
enough. The Fig|2] shows them as distances between the positions of the satellites given by 
the JPL/HORIZONS ephemerides, and the quasiperiodic decompositions of the osculating 
elements, over the 100 years covered by HORIZONS ephemerides. We can see important 
residuals (up to 13,000 km for Epimetheus, while the least Janus-Epimetheus distance is 
about 10,000 km), with a 4-y periodic amplitude, due to the unusual shapes of the orbital 
elements because of the orbital swaps. They induce Gibbs phenomena, i.e. big residuals 
at the swaps. Moreover, there are still residuals far to the swaps, probably due to a lack 
of accuracy in the determination of the short-periodic contributions. These last oscillations 
should be smaller in epicyclic elements. So, we considered this quasiperiodic decomposition 
as not accurate enough to be used in numerical computations. We still are confident in 
the proper modes in TablH this error just means that a large amount of periodic terms is 
required to get accurate enough ephemerides. Moreover, we cannot exclude the presence of 
long-period contributions, for which a polynomial interpolation would be more appropriate 
over this time interval. 

We finally chose to interpolate the cartesian coordinates (i.e. positions and velocities) of 
the two satellit es using the cubic splines interpolation. The reader can find expl anations of 



this method in (iDahlquist fc Bjorckll2008l ). The error associated is of fourth order (ISonneveld 
19691 ). it means that a reduction of the sampling step of the signal by a factor a will reduce 
the interpolation error by a fac tor k, g^. For this, we used the interpolation routines given 



by the GNU Scientific Library (IGalassi et al.l 120091 ). and checked the error in comparing the 



positions given by the interpolation routine with positions given by HORIZONS, at dates 
that were not used to interpolate the ephemerides (at the interpolation points, the error 
should be null) and got an interpolation maximum error of about 850 km with a sampling 
step of 3 hours. We finally used a sampling step of 1 hour to get a maximum error of 9 
km. Contrary to the quasiperiodic decomposition, the error seems to have a quite uniform 
repartition, in particular the dates of the orbital swaps cannot be guessed from the plots of 
errors. We do not show it here because the reader would just see black rectangles. 
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This accuracy of 9 km assumes that the JPL ephemerides are exact over [1950-2050]. 
Their internal error can have at least two causes. The first one is the interpolation error of 
the numerical integration used to compute the ephemerides. This error is considered as lower 
than 25 meters for the DE406 planetary ephemerides, and should be lower for the satellites 
of Saturn, so it is negligible in our study. The other cause of error c omes from the fi t of th e 
dynamical model to the astrometric observations. As explained in iJacobson et al.l (120081 ). 
this fit has been made mostly on Cassini observations since 2004, with a 1 — cr accuracy 
of about 20 km in the downtrack direction. Out of this timespan (i.e. before 2004 and in 
the future) , the accuracy is expe cted to be worse because of the extrapolation. The reader 
can find in iDesmars et al.l (120091 ) a review of the methods used to estimate the accuracy of 
extrapolated ephemerides. 



The rotation 



We here use our interpolation of the JPL HORIZONS ephemerides to get the gravita- 
tional perturbation of Saturn on the rotation of Janus and Epimetheus, with a 1-hour time 
step and so an internal accuracy of 9 km. The influence of the other satellites is accounted 
only as indirect effects, i.e. through perturbations on the orbits of Janus and Epimetheus. 
As for most of the natural satellites of the Solar System, the se two bodies are expected t o 
follow the 3 Cassini Laws, originally enounced for the Moon (jCassinilll693l : IColombdll966l ). 



I.e.: 



1. The Moon rotates uniformly about its polar axis with a rotational period equal to the 
mean sidereal period of its orbit about the Earth. 

2. The inclination of the Moon's equator to the ecliptic is a constant angle (approximately 
1.5°). 

3. The ascending node of the lunar orbit on the ecliptic coincides with the descending 
node of the lunar equator on the ecliptic. This law could also be expressed as: the spin 
axis and the normals to the ecliptic and orbit plane remain coplanar. 

In the case of natural satellites, they can be rephrased this way: the rotation of the satellite is 
synchronous, its angular momentum has a nearly constant inclination on an inertial reference 
plane, and is located in the plane defined by the normal to the orbital plane and to the 
Laplace Plane. The Laplace Plane is the plane normal to the rotation axis of the orbital 
frame, i.e. it is defined with respect to the orbital precessional motion. It has the property to 
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minimize the variations of the orbital inchnations. Using the Laplace Plane as the reference 
frame is of high importance for satellites orbiting far from their parent body, because of the 
Solar perturbation that tends to take the orbital plane away from the equatorial plane of 
the planet. However, for satellites orbiting close to their planet as it is the case here, the 
equatorial plane of Saturn is s o close to the Lapla ce Plane that it c a n be used for des c ribing 
the rotational dynamics (see iDobrovolskid (Il993l ): iNoyelled (120091 ): iTremaine et al.l (120091 ) 
and the Appendix). 



3.1. The model 



The model we used is very similar to the one already used by iHenrardI (j2005al Jbl) for lo 
and Europa and by iNoyelles et al.l (120081 ) for Titan. We consider Janus and Epimetheus as 
rigid triaxial bodies whose matrices of inertia reads 




(7) 



with A< B <C. 



The dynamical model is a 3-degree of freedom one in which 3 references frames are 
considered: 

1. An inertial reference frame (el, €2, 63). We used the one in which the orbital ephemerides 
are given, i.e. mean Saturnian equator and mean equinox for J2000.0 epoch. 

2. A frame {nl, 712,113) bound to the angular momentum of the considered body. 

3. A frame (/i, /2, /s) rigidly linked to the body. 



We first use Andoyer's variables (jAndoyeiill926l : lDepritlll967l ). which are based on two 
linked sets of Euler's angles. The first set {h, K, g) locates the position of the angular 
momentum in the first frame (61,62,63), while the second one, {g,J,l), locates the body 
frame (/i, /2, /a) in the second frame tied to the angular momentum (see Fig. [2D. 

The canonical set of Andoyer's variables consists of the three angular variables l,g,h 
and their conjugated momenta L, G, H defined by the norm G of the angular momentum 
and two of its projections : 
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Date (V) 



Fig. 2. — Distance between the positions given by JPL HORIZONS and those computed from 
our periodic time series of the orbital elements, for Janus (left) and Epimetheus (right). The 
error is due to the difficulty to represent these orbital elements with sinusoidal terms. 




Fig. 3. — The Andoyer variables (reproduced from (iHenrard 
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l L — G cos J 

9 G (8) 

h H = GcosK 

Unfortunately, these variables present two singularities: when J — (i.e., the an- 

— * 

gular momentum is colinear to fz), I and g are undefined, and when K = (i.e., when 
Janus/Epimetheus' principal axis of inertia is perpendicular to its orbital plane), h and g 
are undefined. That is the reason why we shall use the modified Andoyer's variables : 



p=l+g+h P=^ 



r 



-h i?= ^ = P(l -cos/sT) = 2Psin2f ^g) 

^9 = Vising rjg^^cosq 

where n is the body's mean orbital motion , q — —I, and Q — G — L — G{1 — cos J) = 
2G'sin^-|. With these new variables, the singularity on I has been dropped. Using these 
variables has a great mathematical interest, because they are canonical, so they simplify an 
analytical study of the system, as was done in previous works. Our study here is quite purely 
numerical, but wc keep these variables, in order to be consistent with previous studies. We 
later derive other output variables, that arc more relevant from a physical point of view. 

In these variables, the kinetic energy T — ^cu ■ G oi the system reads: 



T 



X 



with 



71+72 _.2 



1-71-72 1-71 + 72 



7i - 72 2 
% 



(10) 



7i 



2C-A-B 
2C 



J, 



C 



(11) 



and 



B-A MU^ 
72 = -1^ = 2(^22-^. (12) 

In these last 3 formulae, u is the instantaneous vector of rotation, M is the mass of Janus or 
Epimetheus, TZ its mean radius, and J2 and C22 the two classical normalized gravitational 
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coefficients related respectively to the oblateness and equatorial ellipticity of the considered 
body. 



The gravitational disturbing potential due to an oblate perturber p reads (IHenrard 



2005cf l: 

= Vpi + Vp2 (13) 

with 

Vpi = -^C^biixl + yl) + ^,{xl - yj)] (14) 

and 



X [li{xl + yl)+l2{xl-yl)l (15) 

where Q is the gravitational constant, Mp the mass of the perturber, its J21 Rp its mean 
radius, dp the distance between the perturber's and Janus (or Epimetheus) ' centers of mass, 
and Xp and yp the two first components of the unit vector pointing to the center of mass of 
the perturber, from the center of mass of the body, in the reference frame {fi, f2, fs)- Vpi 
expresses the perturbation due to a pointmass perturber, while Vp2 represents the perturba- 
tion due t o its J^, assum ing that the body is in the equatorial plane of the perturber. As 



shown in (IHenrardl l2005d ) . it is a good approximation if the sine of the angle between Sat- 
urn's equatorial plane and the orbit is small. Since this quantity is always smaller than 10~^ 
(TabJ23] and [261). we can consider this approximation as satisfying. This assertion assumes 
that the obliquity is very small, what will be checked in this study. 

Usually the orbital ephemerides give us the location of the perturber in the inertial 
frame, so we have to perform 5 rotations to convert the coordinates from the inertial frame 
to (/i, /2, /a). More precisely, if we name (xj, yi, Zi)'^ the unit vector locating the perturber's 
center of mass in the inertial frame, we have 



Xr. 



Vp 



Xi 



R,{-l)R,{~J)R,{-g)R,{-K)R,{-h) y. 



(16) 



Zi 



with 
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cos<p 


— sm 


sin 


cost; 









i?3(0) = I sin0 COS0 I (17) 

and 

/ 1 

= I COS0 —sine; 
\ sin0 COS0 

Finally, the total Hamiltonian of the problem reads: 



H 



X 



7i + 72 



7i 



1-71 -72'' ■ 1 



72 2 

% 



7i + 72 
2 



x[7i(xf^+Z/f^)+72(a;fj- 



(19) 



where the index I] stands for Saturn. We will use this Hamiltonian for a numerical study of 
the rotation. An analytical study can show that the Hamiltonian (fT9!) can be reduced to 



n{u, V, w, U, V, W) = ojJJ + u^V + uj^W + V{u, V, w, U, V, W) 



(20) 



where V represents a perturbation, and the three constants Uu, <^v and uj^ are the periods of 
the free oscillations around the equilibrium defined by the Cassini Laws. This last Hamilto- 
nian is obtained after several canonical transformations, the first one consisting in expressing 
the resonant arguments cr = p — A + vr and p = t + ^ respectively associated with the 1 : 1 
spin-orbit resonance and with the orientation of the angular momentum, A and ^ being the 
orbital variables defined abov e. The complete c alculation is beyond the scope of this paper, 
the reader can find details in jHenrardI boOSal B iNovelles et allboosi ) . 



TabS] gives the values of the physical parameters we used in our numerical integrations. 
They come mainly from spacecraft datas. Pioneer & Voyager for Saturn, and Cassini for 
Janus and Epimetheus. 
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Table 4: Physical parameters of Saturn, Janus and Epimetheus. The gravity parameters J2 
and of Janus and Epimetheus have been derived from the moments of inertia given by 



Tiscareno et al.l (120091 ) . 



Parameters 


Values 


References 


Saturn 




71-^ (mean) 


58232 km 


lAU 2006 fSeidelmann et al. 2007) 


(equatorial) 


60330 km 


Pioneer & Vovaeer fCampbell & Anderson 1989) 




1.6298 X 10-2 


Pioneer & Vovaeer fCampbell & Anderson 1989) 




3.7931272 x 10^ km^.s'^ 


Pioneer & Vovaeer fCamubell & Anderson 1989) 


Janus 




QMj 


0.12660 km3.s-2 


f Jacobson et al. 2008) 


n 


89.5 km 


Cassini fTiscareno et al. 2009) 




0.360 


Cassini fTiscareno et al. 2009) 




0.407 


Cassini fTiscareno et al. 2009) 


C/MU^ 


0.470 


Cassini fTiscareno et al. 2009) 


J2 


8.65 X 10-2 




C22 


1.175 X 10-2 




Epimetheus 




GMe 


0.03513 km3.s-2 


f Jacobson et al. 2008) 


R 


58.1 km 


Cassini fTiscareno et al. 2009) 




0.328 


Cassini fTiscareno et al. 2009) 




0.469 


Cassini fTiscareno et al. 2009) 




0.476 


Cassini fTiscareno et al. 2009) 


J2 


7.75 X 10-2 




C22 


3.525 X 10-2 
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3.2. Numerical integrations 



In order to integrate numerically the system, we first express the coordinates of the 
perturber {x'^,y'^) with the numerical ephemerides and the rotations given in (Eq IT6l) . in 

the body frame (/i,/2,/3)- As explained before, the ephemerides have been obtained in 
interpolating the HORIZONS cartesian coordinates of Janus and Epimetheus with cubic 
splines. This way, we get coordinates depending of the canonical variables. Then we derive 
the equations coming from the Hamiltonian (fT9|) : 



dp 


dH 


dP 


dH 


~(E ~ 


ap' 


dt 


dp ' 


dr 


dH 


dR 


dH 


dt ~ 


~dR' 


dt 


dr ' 


diq _ 


dH 


dr]g 


dH 


dt 


dr]q 


dt 





(21) 



We integrated over 100 years using the Adams-Bashforth-Moulton 10th order predictor- 
corrector integrator. The solutions consist of two parts, the forced one, directly due to the 
perturbation, and the free one, that depends on the initial conditions. The initial conditions 
should be as close as possible to the exact equilibrium, to have low amplitudes of the free 
librations. Several methods exist to reduce the amplitudes of the free librations: 



Bois fc Rambauxl (120071 ) propose to fit the mean initial conditions in order to locate 
the spin-orbit system at its center of libration. 



Peale et al.l (120071 ) add a damping in the equations that reduces the amplitude of the 
free librations. 



Yseboodt fc MargotI (120061 ). in the framework of a numerical integration of the spin 
and of the orbit of Mercury, start from a simple Sun-Mercury system in which the 
equilibrium is very easy to determine analytically, and slowly switch on the planetary 
perturbations in order to create an adiabatic deviation of the equilibrium without 
creation of any free libration. 



NoyellesI (120091 ) use an iterative scheme to remove the free librations from the initia l 
conditions. This last approach is quite similar to the one of iLocatelli fc Giorgilil (120001 ) 
in the framework of a computer-assisted proof of the KAM theorem. 
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The method we used here is inspired from the one of iPeale et al.l (120071 ) . In fact, we 
performed a first numerical integration with a damping, then a reverse integration without 
any damping. This way, we got initial conditions that we use to get the solution of the 
rigid rotation of Janus and Epimetheus without any free rotational energy. We finally use 
these last initial conditions to get a rotational state in which the energy is minimized. We 
also avoid any shift of the equilibrium that a too fast dissipation might induce. The reasons 
why we want to minimize the free component of the solutions is first because this solution 
is expected to have been damped by dissipations, and second because a signal from which 
useless components have been dropped will be more efficiently analyzed. 

Our model of dissipation is just a mathematical one, i.e. without considering the physical 
cause of the dissipation, because it just aimed at finding initial condition and not at studying 
the dissipation itself. We added to every equation a dissipative term alike —a{x — x*) where 
X* is an approximation of the equilibrium of the variable x, and a a positive constant. Our 
only constraint on the value of a is that the dissipation is adiabatic, because a too fast 
dissipation changes the equilibriur a of the system (see for instance the infiuen ce of tides on 
the equilibrium state of Venus in (ICorreia et al.ll2003l : ICorreia fc Laskarll2003l ). respectively 
for the theory and the numerical confirmation). 



3.3. The output variables 

In order to deliver theories of rotation that can be easily compared with observations, 
we chose to express our results in the following variables: 

• Longitudinal librations, 

• Latitudinal librations, 

• Orbital obliquity e (the orientation of the angular momentum of Janus/Epimetheus 
with respect to the normal to the instantaneous orbital plane), 

• Motion of the rotation axis about the pole axis. 

There are at least two ways to define the longitudinal librations. We can for instance 
consider the librations about the exact synchronous rotation, i.e. p— < n > t, usually called 
physical librations. The determination of a "mean" mean motion for Janus and Epimetheus 
is far to be obvious because of the orbital swaps. We can either consider a mean motion 
averaged over several swaps, and get averaged librations that cannot actually be observed 
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at a given date, or consider two different mean motions for each satellite, in assuming that 
the mean motion is constant between two swaps. In order to keep this assumption valid, we 
shall use two intervals of study, necessarily smaller than 4 years, that are far enough from 
the swaps to not be affected by the transitions. We determined the numerical values of the 
mean motions out of the swaps in using the slope of the mean longitudes of each bodies over 
the time spans [1998.5:2001.7] and [2002.5:2005.7] (cf. TabJS]). 



Table 5: Mean motions of Janus and Epimetheus between two orbital swaps, determined in 
using the JPL HORIZONS ephemer ides. These values for the period [2002.5 : 2005.7] are 
close to the mean motions given by iJacobson et al.l (120081 ). fitted on quite the same time 
span. 



Janus 
Epimetheus 



[1998.5 : 2001.7] 



[2002.5 : 2005.7] 



3304.35631 ±5.658 x 10"^ 
3302.78370 ± 1.259 x lO'^ 



3303.67315 ± 5.662 x 10"^ 
3305.24602 ± 1.255 x 10^^ 



Another way to consider the longitudinal librations is to work on the librations about the 
Janus-Saturn (or Epimetheus-Saturn) direction. We will call these librations tidal librations 
because they represent the misalignment of the tida l bulge of the satell i te. T he reader 
can f ind graphical de s cripti ons of these librations in [Murray &: DermottI (Il999l ). Fig. 5. 16 



or m 



Tiscareno et al.l (120091 ). Fig.3, where the physical librations are written as 7 and the 



tidal librations as ip. The difference betwteen the two is the optical lihration (^/^ — 7 = 
2e cosnt), which arise s frorn . Kepler's Third Law. There is a mistake in the Fig. 5.16 
of [Murray fc DermottI (119991 ). that illustrates a dynamically forbidden state (the satellite's 
long axis can never point between the direction towards Saturn an dthe d irection towards 



the em pty focus), the equations being correct. This problem is corrected in lTiscareno et al. 
J2009h . 



The latitudinal librations are the North-South librations of the large axis of the consid- 
ered body in the saturnocentric reference frame that follows the orbital motion of the body. 
They are analogous to the tidal librations that are the East- West librations. In order to 
get to the tidal longitudinal librations and the latitudinal librations, we first should express 
the unit vector /i (i.e. the direction of Janus/Epimetheus' long axis) in the inertial frame 
(61,62,63). From (EqJT6l) and the definitions of the Andoyer modified variables (Eq|9]), we 
get: 



— * 

/i = (cosr(cos(p ± r — /) cos/ — sin(p ± r — /) cos Jsin /) ± 

sin r(cos K{sni{p + r — I) cos / ± cos(p ± r — /) cos J sin /) — sin K sin J sin l))ei 
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+ (— sin r (cos(p + r — I) cos / — sin(p + r — I) cos J sin /) + 

cosr(cosii'(sin(p + r — /) cos / + cos(p + r — /) cos J sin /) — sin K sin J sin /))e2 
+ (sin (sin(p + r — I) cos / + cos{p + r — I) cos J sin /) + cos K sin J sin /)e3. (22) 

The tidal longitudinal librations ip and the latitudinal ones r] are found this way: 



^ = t- h 



(23) 



and 



ri = n-h, (24) 

where n is the unit vector normal to the orbit plane, and t the tangent to the trajectory. We 
get these last two vectors by: 



n 



and 



t 



XXV 

\x X v\ 



n X X 



(25) 



I |n X x| I 

where x is the position vector of the body, and v its velocity. 



(26) 



Finally, the motion of the rotation axis about the pole is derived from the wobble J, it 
is given by the two variables Qi and Q2 defined as: 



Qi = sin J sin/ 1 1 



Jo + 2a 



22 



C 



(27) 



and 



Q2 = sin J cos/ ^1 + 



J2 - 2C. 



22 



C 



(28) 



they are the first two components of the unit vector pointing at the instantaneous North 
Pole of Janus' (respectively Epimetheus') rotation axis, in the body frame of Janus (or 
Epimetheus). These quantities are finally mul tiplied by the polar rad ius of the satellite (76.3 
km for Janus and 53.1 km for Epimetheus (iTiscareno et al.ll2009l )) to get a deviation in 
meters. 



As for the previous study and for the orbital ephemerides, we will give these solutions 
under a semi-analytical (or synthetic) form, thanks to the frequency analysis. 
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3.4. The results 

As said above, there are two parts in the solutions: the free and the forced ones. The 
free solutions depend only on the initial conditions and are assumed to be damped because 
of dissipations (in particular tidal dissipations). It is anyway interesting to study it, at least 
to compute the three periods associated, in case of one of them would be close to a resonance 
with a forced contribution. We detect the free terms thanks to the frequency analysis and 
their values are gathered in TablHl 

Table 6: Periods of the free librations, determined numerically. 





Janus 


Epimetheus 


T 


1.26713 d 


0.74717 d 


T 


2.17884 d 


1.80386 d 


T 


2.09798 d 


5.54234 d 



The forced solutions are gathered in the TabjT] to [TTJ In all the following tables, the 
identification of the periodic contribution has been made in checking the frequencies and 
the phases. In particular, the phases, given at J2000, are very useful to discriminate 2 
contributions with very close frequencies, like wj and we- We here present all the sinusoidal 
terms actually detected by the frequency analysis algorithm, except when we precise in the 
caption that a cut-off has been made. 

Table 7: Tidal longitudinal librations if) of Janus. The phases are given at J2000, the series 
are in sine. 



N 


A 





Wj 


Amplitude 


Phase 


Period 


1 


1 




-1 


1.03° 


116.360° 


0.69735 d 


2 


1 


-1 


-1 


17.54 arcmin 


-61.337° 


0.69752 d 


3 


1 


1 


-1 


—17.52 arcmin 


114.056° 


0.69719 d 


4 


1 


2 


-1 


3.04 arcmin 


112.262° 


0.69702 d 


5 


1 


-2 


-1 


2.99 arcmin 


120.439° 


0.69769 d 


6 


1 


-3 


-1 


— 1.70 arcmin 


123.544° 


0.69785 d 


7 


1 


3 


-1 


1.67 arcmin 


-70.825° 


0.69685 d 



For Janus, the librations are dominated by short-period contributions essentially due 
to the mean anomaly A — w,j for the physical longitudinal librations 7, and to A — for 
the latitudinal ones. The physical cause of these librations is the variations of the distance 
Saturn- Janus because of the eccentricity. We can also notice the other contributions, with 
this form for latitudinal librations: 
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Table 8: Latitudinal librations rj of Janus, in arcseconds. The series are in sine. 



N 


A 







Amplitude 


Phase 


Period 


1 


1 




-1 


-6.18 


-48.717° 


0.69186 d 


2 


1 


1 


-1 


1.79 


-50.972° 


0.69170 d 


3 


1 


-1 


-1 


-1.79 


113.539° 


0.69202 d 


4 


1 


2 


-1 


-0.31 


-53.312° 


0.69153 d 


5 


1 


-2 


-1 


-0.31 


-44.149° 


0.69219 d 



A{ sin(A + i0 - ^j) + sin(A - # - ^j)) = 2A sin(A - ^j) cos(i0) (29) 

where i is an integer. For longitudinal librations, wj should replace Qj in (Eq J2^ . We here 
consider that the amplitudes associated are the same, in fact they have relative differences 
smaller that 10%. Moreover, the phases are consistent with the frequencies, because they are 
contained in the proper modes A, and ^. If the signs of the amplitudes are opposite, the 
right-hand side of Eq. (l29!) is 2A cos(A — ^j) sin(i0), so the resulting wave has the same visual 
aspect. These contributions result in 0.697-d-periodic beatings in 4/i-y-periodic envelopes 
(where i=l,2,3,. . . , so periods of 4, 2, 1.33, . . . years, with decreasing amplitudes associated), 
they are librations due to the orbital swaps. We can see that for Janus, they remain small 
compared to the "classical" librations (i.e. just involving A and the node/pericentre). 

The TabJH] and (TD] present the same variables for Epimetheus. The striking difference 
is that the effects of the orbital swaps dominate the dynamics, whereas the contributions 
A — zue and A — are of smaller amplitude. This difference is due to the amplitude of the 
orbital swap for Epimetheus, that is larger than for Janus because of the mass ratios of the 
two satellites, as we already noticed. 

The results related to the orbital obliquity e are gathered in TabJTTl for the two satellites. 
The phases are not given, but the reader can guess them from Tab{T] and the signs of the 
amplitude (a minus sign means a phase shift of 180°). The most interesting is the first 
line, giving the mean theoretical obliquities, respectively 5.95 and 10.8 arcseconds. Their 
variations, smaller than one arcsec, are too small to be detected from observations. As for 
the librations, the orbital swaps have a bigger effect on Epimetheus than on Janus, but even 
for Epimetheus the obliquity only varies by a few percent. 



As already pointed out by iBouquillon et al.l (120031 ) , there a small motion of the North 



Pole axis about the angular momentum, but also far too small to be detected (^i 1 m for 
the two bodies). The tables associated are in the supplemental material. 

We now give (see Tab JT2|) an estimation of the error induced by the synthetic repre- 
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Table 9: Tidal longitudinal librations of Epimetheus. The series are in sine. 



N 


A 




777 TP 


Amplitude 


Phase 


Period 


1 


1 


— I 


—1 


5 19° 


— 170 621° 


69752 d 


2 


1 


1 


—1 


—5.17° 


4.964° 


69719 d 


3 


1 


2 


-1 


-3.93° 


-177.326° 


0.69702 d 


4 


1 


-2 


-1 


-3.92° 


-168.330° 


0.69769 d 


5 


1 




-1 


-2.58° 


-172.826° 


0.69735 d 


6 


1 


3 


-1 


-1.09° 


0.215° 


0.69685 d 


7 


1 


-3 


-1 


1.07° 


-165.862° 


0.69785 d 


8 


1 


-4 


-1 


28.14 arcmin 


15.891° 


0.69802 d 


9 


1 


4 


-1 


27.50 arcmin 


-1.553° 


0.69669 d 


10 


1 


5 


-1 


20.03 arcmin 


175.657° 


0.69652 d 


11 


1 


-5 


-1 


— 19.56 arcmin 


18.683° 


0.69819 d 



Table 10: Latitudinal librations rj of Epimetheus, in arcseconds. The series are in sine. 



N 


A 







Amplitude 


Phase 


Period 


1 


1 


1 


-1 


5.75 


-89.666° 


0.69170 d 


2 


1 


-1 


-1 


-5.75 


94.847° 


0.69202 d 


3 


1 


2 


-1 


4.49 


88.082° 


0.69153 d 


4 


1 


-2 


-1 


4.48 


97.010° 


0.69219 d 


5 


1 




-1 


2.68 


92.593° 


0.69186 d 


6 


1 


3 


-1 


1.30 


-94.153° 


0.69137 d 


7 


1 


-3 


-1 


-1.30 


99.340° 


0.69235 d 



sentation of the output variables. This is the maximum error over the time interval of the 
study (i.e. [1950:2050]), that we compare with the maximum amplitude of the variable. In 
comparing this maximum error with the amplitudes of the sinusoidal terms given in Tab 17] to 
[TTl it seems to be too high. In fact, this maximum error is due to the Gibbs phenomenon at 
the orbital swaps that comphcates the convergence of the periodic series, as for the orbital 
motion (FigjS]). Moreover, there is also an error far from the swaps, probably due to a lack 
of accuracy in the determination of the short periods. The FigJH illustrates this problem for 
the longitudinal librations. It also shows that the error is smaller at the center of the interval 
1 arcmin) than at the edges (~ 2 arcmin), probably because of long-period terms that 
the frequency analyses did not identify. 

As explained earlier, we also computed the physical longitudinal librations (FigJS]) on 
the time intervals [1998.5 : 2001.7] and [2002.5 : 2005.7], i.e. in excluding the orbital swaps. 
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Table 11: Orbital obliquity e of the two bodies, in arcseconds. The contributions for which 
no amplitude is given have not been detected by the frequency analysis. It was required to 
remove a slope of —5.799 x 10~^ arcsec/y, probably the signature of a long-period contri- 
bution, zu stands for zuj for Janus, zue for Epimetheus. Contributions associated with the 
nodes ^ have been actually detected by the frequency analysis, but with smaller amplitudes 
(< 0.01 arcsec for Janus, and < 0.3 arcsec for Epimetheus). The series are in cosine. 



A 




to 


Amplitude 
Janus 


Amplitude 
Epimetheus 


Period 








5.94 


10.83 


oo 


1 




-1 


0.02 


-0.43 


0.69735 d 


1 


-1 


-1 


0.01 


0.87 


0.69752 d 


1 


1 


-1 


-0.01 


-0.87 


0.69719 d 


1 


2 


-1 




-0.66 


0.69702 d 


1 


-2 


-1 




-0.66 


0.69769 d 



in considering two different values of the mean orbital motions (Tab 15]). The most striking 
is here the presence of a long-period contribution (probably 8 years), showing the influence 
of the orbital swap on the rotation. We can also see some short-periodic variations of the 
amplitudes. Unfortunately, a frequency analysis over such a short time span cannot split 
contributions alike A ± — from A — cc, whereas it was possible for the tidal librations 
over a wider interval of study (TabJ7]and[9]). 

In fact, a frequency analysis would only give a kind of averaged short-period libration 
(TabJT3!). The variations of the amplitudes associated can be guessed in reading the plots. 

3.5. Analytical validations 

Many analytical studies exist on the synchronous rotation. We here propose to use them 
to validate our numerical results. 

• Mean obliquity 

As explained in the Appendix, the Third Cassini Law implies that the normal to 
the reference plane (that should be an acceptable Laplace Plane), the normal to the 
orbital plane and the angular momentum of Janus / Epimetheus are coplanar. As a 
consequence, 4 equilibria known as Cassini States are possible, in a simplified model. 
Another analytical formula for the Cassini State 1 (i.e. the most probable one) is given 
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Table 12: Maximum amplitudes and errors of the quasiperiodic representation of the output 
variables of the rotation, estimated from the plots. The error is mostly due to the Gibbs 
phenomenon at the orbital swaps, see also FigJH 



Variable 


Amplitude Error 




Tidal longitudinal librations 


Janus 


^1.1° ^5 arcmin 


Epimetheus 


^ 10° ^ 1.5° 




Latitudinal librations 


Janus 


~ 6.5 arcsec ~ 0.9 arcsec 


Epimetheus 


^11 arcsec ~ 3 arcsec 




Obliquity e 


Janus 


^6.15 arcsec ^0.2 arcsec 


Epimetheus 


^ 14 arcsec ~ 1 arcsec 



Table 13: Longitudinal physical librations 7. 





[1998.5 


: 2001.7] 


[2002.5 


: 2005.7] 




Frequency 


Amplitude 


Frequency 


Amplitude 




rad/y 




rad/y 




Janus 


3291.2598 


0.337 ±0.03° 


3290.5893 


0.338 ±0.03° 


Epimetheus 


3289.7159 


8.615 ±0.9° 


3292.1325 


8.588 ±0.9° 



by iHenrard fc SchwanenI (120041 ). In that paper, the equilibrium obliquity K* related 
to the orbital plane is given by: 



with 



and 



^ _Jl±^i (30) 



di = -^2—^ (31) 



62 = -SC22 ^ , (32) 



what gives straightforwardly, using e = K* — 1: 



1950 1960 1970 1980 1990 2000 2010 2020 2030 2040 2050 

Date (V) 



1950 1960 1970 1980 1990 2000 2010 2020 2030 2040 2050 

Date (vl 




Fig. 4. — Tidal longitudinal librations ip of Janus (up) and Epimetheus (down). The panels 
(b) and (d) show the residuals with the quasi-periodic representations, on which we can see 
the signature of a 4-y periodic error, probably due to the Gibbs phenomenon at the orbital 
swaps. 



^ ~ 1 _L 3 n J2+2C22 • (^^) 
2 ^ C/M-R? 

This last formula is very similar to the one given in the Appendix for the Cassini State 
1 in assuming that sin / I and cos / ~ 1, which holds for / << 1. 

The table [TH gathers the locations of the theoretical Cassini States and recalls the 
mean obliquities that we get from our numerical simulations. We can notice that our 
results are close to the Cassini State 1. We get for Janus significant differences between 
our values and the analytical ones, while the agreement is very good for Epimetheus. 
These values have been computed in neglecting the 1 : 1 orbital resonance, and in 
assuming that the orbits of the satellites were circular and uniformly precessing. A 
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1998.5 1999 1999.5 2000 2000.5 2001 2001.5 



2002.5 2003 2003.5 2004 2004.5 2005 2005.5 



□a1e (y) 



□a1e (y) 





2000 2000.5 

□a1e (y) 



2004 2004.5 

□a1e (y) 



(c) (d) 

Fig. 5. — Physical longitudinal librations 7 of Janus (up) and Epimetheus (down). These li- 
brations have been obtained in removing a fitted mean spin from the orientation of the longer 
axis on the inertial reference frame, respectively 3304.356301 rad/y (a) and 3303.673150 
rad/y (b) for Janus, and 3302.7839016 (c) and 3305.24602 (d) for Epimetheus. 



possible explanation of the discrepancy observed for Janus could be the influence of a 
long period that required to fit a slope on the obliquity of Janus (TabJTTl). 

Fundamental frequencies of the free librations 

We can use previous studies to esti mate analytically the three proper frequencies of 
the librations about the equilibrium. iHenrard fc SchwanenI (120041 ) detailed an extensive 
derivation of the Hamiltonian (fT9|) for the satellites in 1 : 1 spin-orbit resonance, in 
assuming that the orbit of the satellite was circular and uniformly precessing, but 
in considering the three degrees of freedom of the system. After determination of the 
equilibrium and several canonical transformations, they got the following Hamiltonian: 



M{u, V, w, U, V, W) = ujJJ + uj^V + uJwW, 



(34) 
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Table 14: Theoretical Cassini States for Janus and Epimetheus, compared with our numerical 
simulations. The calculations use the equations given in the Appendix. 





Janus 


Epimetheus 




1161.73 


1542.84 




-89.099 


-118.315 




-1.030 


-1.051 


Cassini State 1 


6.719 arcsec 


10.834 arcsec 


Cassini State 2 


-89.357° 


-89.516° 


Cassini State 3 


179.998° 


179.997° 


Cassini State 4 


89.357° 


89.516° 


Numerical simulations 


5.945 arcsec 


10.829 arcsec 



in which the constants a;„, and oj^ are the free frequencies of the small librations 
about th e equilibrium. This analytical method has been successfully used in (IHenrard 
2005al lbl: iNm^elles et all 120081 ) 

In a more simplified model (i.e. no wobble and a negligible obliquity), we can con- 
sider that the lib r ation angle 7 about the synchronous rotation is ruled by (e.g. 
Goldreich fc Peale JigS)): 



+ -n^{B - A)H{l,e) sin 27 = 0, 



(35) 



with -ff(l, e) = 1 — 5e^/2 + 13e^/16 + O(e^). Since 7 and the eccentricity e are small, 
the equation fl35l) becomes: 



with 



(36) 



B-A 



3C. 



The equation (136|) is known to rule a pendulum swinging at the frequency cu^, i.e. the 
frequency of the free longitudinal librations. We propose in the Appendix a derivation 
of the free librations of the obliquity. 

The TabJT^ gives a comparison between our numerical values of the periods of the 
free librations and values due to formulae obtained by analytical studies. The column 
"Analytical 2" uses the Eq. (lC7p that uses the value of obliquity coming from our 
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Table 15: Comparison between the periods (in days) of the free hbrations obtained numeri- 
ca lly and in usin^ analyt i cal fo rmulae. The column "HS 2004" refers to the analytical works 
of iHenrard fc SchwanenI (120041 ). while "Analytical 1" refers to Eg. (1371) for T„ = I-r jiOu and 
to Eq.dCe]) for T^. "Analytical 2" refers to to Eq.dCT]) for T^. 





Numerical 


HS 2004 


Analytical 1 


Analytical 2 


Janus 




T 


1.26713 


1.26053 


1.26814 




T 


2.17884 


1.98573 


1.99787 


1.76767 


T 


2.09798 


1.74391 






Epimetheus 




T 


0.74717 


0.73255 


0.73697 




T 


1.80386 


1.49107 


1.50015 


1.49946 


T 


5.54234 


4.71920 







simulations. We can see a good agreement for the period of the free librations in 
longitude T„, but a significant discrepancy for the two other periods, that the use of the 
real obliquity does not reduce. This discrepancy is probably due to the approximation 
on the dynamics of the satellites, in particular the effects of the 1 : 1 orbital resonance. 

Longitudinal librations 

We have defined two k inds of longitudinal l i bratio ns: the physical ones 7, and the tidal 
ones We have (e.g. iMurray &: Dermott Jl999h ): 



2e 



7 



sin(nt + ip), 



(3^ 



1 - (n/cj„)2 

where cp is a constant phase depending on the time origin. This formula assumes that 
the orbit is close to be keplerian and circular, in particular its eccentricity is assumed 
to be small and constant. The tidal librations tp, representing the misalignment of 
the tidal bulge of the satellite, is obtained in considering the optical libration = 
2esinnt = 7 — 'i/^, due to the var iations of the velocity of the satellite on its orbit. We 
then find (ITiscareno et al.l 120091 ): 



-2e 



1 - {uju/ny 



sin(nt + (p). 



(39) 



The TabJT6] gives a comparison between the amplitudes of the librations that we de- 
tected numerically, and the ones given by the analytical formulae (155]) and (15^ . in 
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Table 16: Comparison between the numerical and analytical determinations of the librations 
of Janus and Epimetheus. The column "Numerical" gathers the results given by the TabJT2] 
and [131 while the "Analytical" ones use the formulae fl38l) and fl39l) with two different values of 
the proper frequency Uu'- the numerical one in column "1", and the one given by the formula 
( 157]) in the last column. These calculations have been made using the eccentricities at J2000.0 
that we determined in fitting a slope, i.e. 7.30539 x 10~^ for Janus and 1.01621 x 10~^ for 
Epimetheus. 





Numerical 


Analytical 1 


Analytical 2 


Janus 




7 


^ -0.34 ±0.03° 


-0.35959° 


-0.35877° 


V. 


^ -i.r 


-1.19673° 


-1.19591° 


Epimetheus 




7 


-8.6 ±0.9° 


-7.41096° 


-9.26011° 




^ -10° 


-8.57545° 


-10.42460° 



using the proper period detected numerically (column "1") and analytically (col- 
umn "2"). For Janus, we can see a quite good agreement between the numerical and 
analytical values, even if the discrepancy is significant. However, an agreement is not 
so obvious for Epimetheus, even between the analytical values. The reason is proba- 
bly a high sensitivity on the gra vity field parani e ters ( J?, C22 and C/MTZ"^) probably 
because, as already indicated by lTiscareno et al.l (120091 ) . Epimetheus seems to be close 
to the 1 : 1 secondary resonance between its spin and its free longitudinal librations. 



Discussion 



Recently, iTiscareno et al.l (120091 ) observed the rotation of Janus and Epimetheus thanks 
to the Cassini spacecraft. We here propose to use the theory of the rotation we have elab- 
orated to try to explain the observations, and to get planetological consequences that could 
be observed. 



4.1. Uncertainty of the solutions 



All our calc ulations assume that t he input parameters are known with a very high 
accuracy. In fact, ITiscareno et al.l (120091 ) derived them from fits to observations, giving most 
probable values and uncertainties (see Tab ITT]) . A rigorous way to study the influence of 
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the uncertainties of the input parameters on the rotation would be to run the numerical 
integrations with some variations of these parameters within the error bars. This method 
requires lots of CPUs and would give uncertainties on useless values, like amplitudes of some 
periodic contributions that are far too small to be observed. That is the reason why we 
propose to use the analytical formulae to derive uncertainties on the longitudinal librations 
and the equilibrium obliquities. 

For this, two parameters need to be known for each body: {B — A)/C for the amplitude 
of the longitudinal librations (EqJSSl and 15^. and (C — A)/C (EqJB2D for the equilibrium 
obliquity. The uncertainties on {B — A)/C are given by iTiscareno et al.l (120091 ) and so can 
be used directly, while the ones on {C — A)/C should be guessed from the uncertainties on 
the dimensions of the bodies. In assuming that they are homogeneous ellipsoids, we have 



what yields 



A 



M 



{b' + c'),B 



M 



M 



(40) 



C-A 
C 



+ 62- 



(41) 



Unfortunately, this last formula does not work for Epimetheus because, as lTiscareno et al. 



( 120091 ) noticed, it in fact deviates s ignifi cantly from an ellipsoid. Moreover, the moments 
of inertia given bv ITiscareno et al.l (20091) are taken directly from a numerical shape model 



derived from the data ( ThomaJ 



1993 



Simonelli et al.lll993l : lThomas et al.lll998l ). and are not 



directly related to the best-fit ellipsoid. We also had troubles in using it, getting obliquities 
far bigger than the mean obliquity (10.8 arcsec) we give in Tab ITT] or the maximum obliquity 
14 arcsec (TabEl). So, we used it only for Janus and got (C - A)/C = 0.234t[|;[]i. The 
Tab ITS] gathers our results. 



Tiscareno et al.l (120091 ) predicted the physical librations 7 as —0.33 ± 0.06° for Janus 



and —8.91^5^4 ° for Epimetheus, from shape-derived moments of inertia that assume constant 
density. We have good agreements, the differences being likely to be due to different con- 
sidered values for the eccentricities. As already said, the longitudinal motion of Epimetheus 



Table 17: Parameters derived from Cassini observations (Tiscareno et al. 2009). 




{B -A)/C a (km) b (km) c (km) 


Janus 

Epimetheus 


0.100 ±0.012 101.5 ± 1.9 92.5 ± 1.2 76.3 ±1.2 
0.296l[J:[!^? 64.9 ±2.0 57.0 ± 3.7 53.1 ± 0.7 
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Table 18: Uncertainty on the longitudin al librations and on th e obliquity, deduced from the 
uncertainty on the parameters given by iTiscareno et al.l (120091 ) . 



7 [1 

Mean obliquity (arcsec) 



Janus Epimetheus 



-0.36 ±0.06 
-1.2 ±0.06 
6.72 ±0.82 



-9.2 



-4.3 
-10.8 



-in 4+4-4 



is very sensitive to the input parameters because of the proximity of the 1 : 1 secondary 
resonance between the spin and the free longitudinal librations. This way, when the input 
parameter [B — A)/C gets closer to the critical value 1/3, the amplitude of the libration 
increases dramatically (Figl6]). We have a similar phenomenon for t he longitudinal lib rations 
of Mercury, close to a resonance with the mean motion of Jupiter (jPeale et al.ll2009l ). 




0.265 0.27 0.275 0.28 0.285 0.29 0.295 0.3 0.305 0.31 0.315 

(B-A)/C 



Fig. 6. — Dependency of the longitudinal librations of Epimetheus on the moments of inertia. 
We can notice a high sensitivity due to the proximity of a resonance. 
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4.2. Comparison with the observations 



The main results of iTiscareno et al.l (120091 ) on the rotation of Janus and Epimetheus are 
detections of the physical librations of Janus and Epimetheus (cf. Tab JT9|) . and a significant 
offset of Janus' sub-Saturn point from the axis of the minimum moment of inertia. 



Table 19: Comparison between the observations of ITiscareno et al.l (120091 ) and our results. 
The column "our simulations" refers to the results of the numerical integration (TabJT3]and 
TTl) while mean values refers to the study of the uncertainties (Tab JTSj) . 





Tiscareno et al. (2009) 


Our simulations 


Mean values 


Janus 




7 


-0.3 ±0.9° 


-0.34 ±0.03° 


-0.36 ±0.06° 


Obliquity e 




5.95 arcsec 


6.72 ±0.82 arcsec 


Longitudinal offset 


5.2 ± r 







Latitudinal offset 


2.3 ± 1° 







Epimetheus 




7 


-5.9 ± 1.2° 


-8.6 ±0.9° 


_q q+4.9 o 


Obliquity e 




10.83 arcsec 




Longitudinal offset 


< 1° 







Latitudinal offset 


< 1° 








Our study seems to agree with the observed value of the physical longitudinal librations 
7 for Janus, even if we should keep in mind that the error bars of the observations are thrice 
bigger than the most probable value, so we cannot be certain that this libration has actually 
been observed. For Epimetheus, there is a significant discrepancy on the mean values of 7, 
but with an overlap due to the sensitivity on the input parameters. An interesting point is 
the error bar coming from our numerical simulations, that assumed the input parameters to 
be exactly known. This error bar is due to the variations of the amplitude that we read on 
th e plots (^FigJSl). The va riation of 0.9° could partly explain the error bar of 1.2° observed 
by ITiscareno et al.l (120091 ). The expected obliquities are probably too small to be detected. 



Our model does not explain the offset of Janus' sub-Saturn point from the axis of the 
moment of inertia. If we suppose it could be due to the orbital swap or to the latitudinal 
librations, it should anyway be lower than 40 arcmin in longitude (FigJSj) and that 6.5 arcsec 
in latitude (Tab JT^ . One way to get it might be to introduce non-diagonal terms in the 
matrix of inertia (Eq|7]). This would physically mean that the geometrical principal axes 
differ significantly from the gravitational ones, so that Janus is not an homogenous body, or 
that its equipotential surface differs significantly from an ellipsoid. We leave this for future 
work. 
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5. Conclusion 

We have elaborated in this paper a 3-degree of freedom theory of the rotatiori of Janus 



and Epimetheus, based on the physical parameters given by lTiscareno et al.l (120091 ) . derived 
from Cassini observations. This study used an accurate representation of the orbital motions 
of these bodies, so as to evaluate the rotational consequences of the orbital swaps. 

Our numerical integrations show that these swaps induce significant variations of the 
short-period longitudinal librations, and thus could explain the error bars on the observa- 
tions, at least for Epimetheus. On the contrary, this model does not explain the mean 
orientation of Janus' axis of minimum moment of inertia. Moreover, we give an estimation 
of some undetected aspects of the rotation, like the latitudinal librations (~ 6.5 arcsec for 
Janus and ^11 arcsec for Epimetheus) and the mean obliqui ties (respectively 6 and 11 
arcsec). Finally we highlight, as already iTiscareno et al. J2009h did, the sensitivity of the 



longitudinal motion of Epimetheus on the input parameters, because of the proximity of a 
resonance. 

It will be probably possible, in a next future, to get a better accuracy on the observations, 
particularly thanks to next fiy-bys. We hope that it will open the door to an inversion of 
the theory of the rotation to get several parameters of the gravity field of these bodies, and 
so information on their interior. 
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A. Approximation of the mean Laplace Plane 

The Lagrange equation ruling the precessional motion of Janus/Epimetheus is, at the 
first approximation (i.e. in neglecting the effects of eccentricities and inclinations) 

+ + ^^^^ 

i 

with 
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and 



2 



cut 



3 f^r^ 

2 \ a 



1 M,; 



W,; = 



4M, 



a,; 



(A2) 
(A3) 



1 Mf. 







4Mi 



[see e.g. 



(— 



^3/2 



i0 



(A4) 



Champenois fc Viennd (119991 )). where n is Janus/Epimetheus' mean motion, 



the semimajor axis of sateUite i (1 to 8 standing respectively for Mimas to lapetus, the 
other ones being neglected), Mj its mass, a© and Mq the same quantities for the Sun, 71^ 

is a cl assical Laplace coeffici e nt. T he Laplace 



the equatorial radius of Saturn, and b^^^^i^) 
coefficients are computed with the formula (see e.g. iBrouwer fc Clemencd ( 119601 )): 



.(1) 
"3/2 



x] 



We have (lYseboodt fc MargotJ 120061 ) 



cos^ 



vr Jo (1 - 2xcos6' + x2)3/2 



d9. 



(A5) 



(A6) 



and 



= : : , (A7) 

where pk = sin^^sinifc and qk = cos sin for the satellites 1 to 8, the Sun and Saturn, 
ik standing for the inclination of the satellite k, and for its ascending node, po and go 
indicate the location of the Laplace Plane. The numerical values of u are gathered in TabJ20l 

A rigorous determination of the instantaneous Laplace Plane would require to consider 
the time variations of the quantities p and q of each perturber, to get po{t) and ^0(^)5 that 
could be finally averaged to get a mean Laplace Plane, that would be an optimized inertial 
reference plane. Since we just want to show that the mean Laplace Place is very close to the 
Saturnian equator at J2000, we will give an upper bound for the inclination of the Laplace 
Plane io, in considering that 
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Table 20: Determination of the Laplace Plane of Janus and Epimetheus. The gravity field 
data and mean semi-major axes come from JPL HORIZONS at J2000, except QMq that 
come from the recommendations of lERS. The u have been computed thanks to Eq lA3l 
and IA4[ The semimajor axes are given with Saturn as the central body. The inclinations 
with resp ect to Saturn's equatorial plane are estimated f rom the main terms of the series of 
TASS1.7 (jVienne fc Duriezlll995l : iDuriez fc Viennd 119971 ). this rough estimation being here 
accurate enough because we just want to show that the inclination of the Laplace Plane is 
small. These values are easier to use than HORIZONS' because HORIZONS give them with 
respect of the Laplace Plane of the considered body. 





GM 


a 


UJ 




sinifc 






km? .s~'^ 


km 










Mimas 


2.5026 


185539 


-9.15 X 10- 


-7 


2.4 X 10- 


-2 


Enceladus 


7.2027 


238037 


-4.01 X 10- 


-7 


2.5 X 10- 


-4 


Tethys 


41.2067 


294672 


-8.03 X 10- 


-7 


1.6 X 10- 


-2 


Dione 


73.1146 


377415 


-5.27 X 10- 


-7 


3.3 X 10- 


-4 


Rhea 


153.9426 


527068 


-3.44 X 10- 


-7 


5.9 X 10- 


-3 


Titan 


8978.1382 


1221865 


-1.41 X 10- 


-6 


1.1 X 10- 


-2 


Hyperion 


0.3727 


1500934 


-3.12 X 10- 


-11 


1.2 X 10- 


-2 


lapetus 


120.5038 


3560851 


-7.44 X 10- 


-10 


0.26 




Sun 


1.32712442076 x 10" 


1.4266641409 x 10^ 


-1.27 X 10- 


-8 


0.47 




Saturn 


3.7931208 X 10^ 





-12.71 










sinio < 



sin iiUJi + sin IqUJq + sin i\^0J^ 

Y,. UJi + UjQ+UJ^ 



(At 



and we get numerically io < 4.63 x 10 rad, i.e. ~ 1 milli-arcsec. 



B. The Cassini states 



Colombo! (119661 ) showed that the spin axis of a rotating body has 2 or 4 equilibria named 
Cassin i states. We give here an analytical study of the location of these states, widely inspired 

but rewritten our way. The reader can find alternative 



(!l987h : iD'Hoedt et al.l f!2006h ) 



Peak 


( 


1969 


); 


Beletskii 


( 


1972): 



We here consider a reference plane of normal k, and we assume that the orbital plane 



- 38 - 



has a constant inclination / and precesses at a uniform rate Under these assumptions, 
the reference plane can be considered as a Laplace Plane. We name n the normal to the orbit 
and s the unit vector colinear to the angular momentum of the considered body. Finally, 
we call e the obliquity of the body, it is the angle between s and n. The Third Cassini Law 
tells us that the vectors k, n and s are coplanar, so we say that the angle between k and s is 



T + e. This convention on the orientation of e is consistent with our model, IWard fc Hamilton 



(120041 ) chose another one. 

The equation of motion ruling s reads: 



^ = a(s ■ n)(s X n) + X fc), (Bl) 
at 



where a is a precessional constant that can be written (IWardlll975l ): 



^ 3 {C - A)r? ^ lr^J2 + 20n , . 

" 2 Co; 2uj C/MR^ ' ^ ^ 

where n is the orbital mean motion and uj is the spin velocity of the body. For most of the 
natural satellites like Janus and Epimetheus, the synchronous rotation induces n = u, and 
we get 

3 J2 ~\~ 2C22 /T-l \ 

At the equilibrium, we have ^ = in the reference frame defined by the vector n, k 
and their product. In fact, n and s precess synchronously. The projection of the equation 
fIBll) on the direction normal to the plane (s, k) yields: 

^sin2e + sin(/ + e) = 0. (B4) 
2^2 

It is well known that, for most of the natural satellites and Mercury, the slow precessional 
motion induces the existence of four equilibria that are close respectively to 0, |, vr, — |. 
These equilibria are known as Cassini States 1, 4, 3 and 2. With a more rapid precessional 
motion, three of these equilibria get closer (except Cassini State 3), and two of them can 
vanish as it is the case for the Moon. 



As can be seen on FiglTJ the value of a;/f2 is critical. It can be shown (IWard fc Hamilton 



20041 ) that the 4 Cassini States exist if \a/^\ > (sin^^^ J + cos^/^ 1)^/^ ^ 1. For the "classical'' 
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Obliauitv e (deal 
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Obliauilv e fdea) 



Fig. 7. — Locations of the Cassini States from the equation (1B4I) . respectively with a/Cl = 
—100 (left), and a/Cl = —1-3 (right), for an orbital inclination / = 0.01 rad. The y axis 
is the quantity a;/(2f2) sin2e + sin(/ + e), the location of the Cassini States being its roots. 
The 4 Cassini States are respectively, from left to right : 1, 4, 3 and 2. We can see that for 
a high value of the ratio |Q;/f2| the four equilibria are close to multiples of 7r/2, while the 
states 1, 4 and 2 get closer when this ratio tends to 1. 

case of the natural satellites where \a/^\ 3> 1, it is straightforward to develop the equation 
(]B4p around 0, |, vr, — | and we get the following formulae: 



sin / 



e ^ — 



r I cos / 



2 

71 + 



sin/-a/^' 
sin / 



cos I 



- + 

2 sin/+Q/^ 



(Cassini State 1) 
(Cassini State 2) 
(Cassini State 3) 
(Cassini State 4) 



C. Free librations of the obliquity 

The aim of this section is to estimate the period of the small free libration of the 
obliquity about the exact Cas sini State 1. We start from the following Hamiltonian (jWard 



1975 



Ward fc Hamilton! I2OO4J ) . representing the behavior of the obliquity, that is in fact a 



part of the Hamiltonian of the system: 



n 



a 



(CI) 



where a, s, n and k are defined as above. We then find: 



a 



H = cos^ e - ^ cos(/ + e) 



(C2) 
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and, at the second order in e: 

H^e^sin/ + -(a + ^cos/). (C3) 

In order to study the small variations of the obliquity e, we set e = eo + ei, where eo is 
the constant obliquity at the equilibrium, and |ei| ^ |eo|. We can now write, in dropping 
the constant terms: 



^ = ei(eo(a + ^cos/) + ^sin/) + ^(a + f^cos/) 



e2 



i(a + ^cos/) (C4) 

because the definition of the Cassini State 1 gives eo(a + ^cosJ) + ^sinJ = 0. We have 
now a quadratic function of ei. 

ei is not really a canonical variable. It results from the classical polar transformation 
of canonical variables that the Hamiltonian flC4p can be written as 



niv,V)=iu^V (C5) 

where v and V are canonical angle-action variables, V = ^ cjt, = a + Q^cosI is the 
frequency of the free oscillations. So, the period of the free librations of the obliquity is, in 
this simplified model: 

n = (C6) 

a + cos I 

and can also be written as 

r„ = (C7) 



D. Supplemental Material 
D.l. Orbital motion 

Here are the frequency analyses of the orbital elements of Janus and Epimetheus. 
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Table 21: Mean longitude of Janus, after removal of the proper mode A. The series are in 
sine. 



N 


cp 


^' Amplitude (rad) 


Period 


1 


1 


0.55346 


8.00687 y 


2 


3 


-0.06351 


2.66892 y 


3 


5 


0.02017 


1.60135 y 


4 




2 0.00798 


68.12 y 



Table 22: z = k + \J —\h = e exp{\^ —Iw) for Janus. The series are in complex exponentials. 



N 


A 





zuj Amplitude (xlO^) 


Period 


1 






1 6.64720 


175.30788 d 


2 


1 




3.59821 


0.69459 d 


3 


1 


1 


1.03092 


0.69442 d 


4 


1 


-1 


1.03032 


0.69475 d 


5 


1 


2 


0.17933 


0.69426 d 


6 


1 


-2 


0.17920 


0.69492 d 


7 


1 


3 


0.09885 


0.69409 d 



Table 23: ( — q + -s/^p = sin exp(-\/^f2) for Janus. The series are in complex 
exponentials. 



N 


cu Amplitude (xlO^) 


Period 


1 


1 - 1.43593 


-176.00968 d 


2 


- 1 0.00797 


136.16498 y 
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Table 24: Mean longitude of Epimetheus, after removal of the proper mode A — tt. The series 
are in sine. 



N 


o 


Aiiii)litu(Io (rad) 


Period 


1 


1 


-1.99423 


8.00676 y 


2 


3 


0.22886 


2.66892 y 


3 


5 


-0.07271 


1.60136 y 


4 


7 


0.03345 


1.14382 y 


5 


9 


-0.01828 


324.94152 d 



Table 25: z — k + \^^h — eexp(-\/^ro) for Epimetheus. The series are in complex 
exponentials. 



N 


A 




Amplitude (xlO''^) 


Period 


1 




1 


9.82435 


175.28140 d 


2 


1 


1 


-2.07104 


0.69442 d 


3 


1 


-1 


-2.07007 


0.69475 d 


4 


1 


2 


-1.59186 


0.69426 d 


5 


1 


-2 


-1.59002 


0.69492 d 


6 


1 




-0.99815 


0.69459 d 


7 


1 


3 


-0.45071 


0.69409 d 


8 


1 


-3 


-0.44971 


0.69508 d 


9 


1 


4 


-0.18343 


0.69393 d 



D.2. Motion of the North Pole 
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Table 26: C = + = sin ^| j exp(V— 1^) for Epimetheus. The series are in complex 

exponentials. 



N 





uj Amplitude (xlO^) 


Period 


1 


1 


3.07916 


-175.99124 d 


2 


-1 1 


0.02859 


-166.00121 d 


3 


1 1 


0.02841 


-187.26047 d 


4 




1 0.00797 


136.15656 y 



Table 27: Motion of the North Pole of Janus' rotation axis about the North Pole of Janus, 
in meters. The variable represented here is Qi + ^/^^Q2 = sin J( sin /(I + (J2 + 2C22)/C) + 
a/— Tcos/(l + (J2 — 2C22)/C)) (see EqJ271 and . The series are in complex exponentials, 
and the cut-off is at 10 cm. 



N 


A 







Amplitude 


Phase 


Period 


1 


1 




-1 


1.16 


-138.717° 


0.69186 d 


2 


-1 




1 


0.79 


-41.283° 


-0.69186 d 


3 


1 


1 


-1 


0.32 


-140.972° 


0.69170 d 


4 


1 


-1 


-1 


0.32 


43.539° 


0.69202 d 


5 


-1 


-1 


1 


0.23 


-39.029° 


-0.69170 d 


6 


-1 


1 


1 


0.23 


136.461° 


-0.69202 d 



Table 28: Motion of the North Pole of Epimetheus' rotation axis about the North Pole of 
Epimetheus, in meters. The series are in complex exponentials, and the cut-off is at 20 cm. 



N 


A 





We 




Amplitude 


Phase 


Period 


1 


1 


1 




-1 


0.72 


-179.666° 


0.69170 d 


2 


1 


-1 




-1 


0.72 


4.847° 


0.69202 d 


3 


-1 


-1 




1 


0.70 


-0.335° 


-0.69170 d 


4 


-1 


1 




1 


0.70 


175.154° 


-0.69202 d 


5 






-1 


1 


0.59 


4.581° 


-87.81766 d 


6 


1 


2 




-1 


0.56 


-1.918° 


0.69153 d 


7 


1 


-2 




-1 


0.56 


7.100° 


0.69219 d 


8 


-1 


-2 




1 


0.55 


-178.082° 


-0.69153 d 


9 


-1 


2 




1 


0.55 


172.900° 


-0.69219 d 


10 






1 


-1 


0.43 


175.419° 


87.17665 d 


11 


1 






-1 


0.34 


2.593° 


0.69186 d 


12 


-1 






1 


0.33 


177.407° 


-0.69186 d 
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